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DEUTERON DISTRIBUTION IN NUCLEI AND IN CORRELATED 

NUCLEAR MATTER 
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We compute the distribution of quasideuterons in doubly closed shell nuclei and infinite 
correlated nuclear matter. The ground states of 16 O and 40 Ca are described in Is coupling 
using a realistic hamiltonian including the Argonne v' s and the Urbana IX models of two- 
and three— nucleon potentials, respectively. The nuclear wave function contains central and 
tensor correlations, and correlated basis functions theory is used to evaluate the distribution 
of neutron-proton pairs, having the deuteron quantum numbers, as a function of their total 
momentum. By computing the number of deuteron-like pairs we are able to extract the 
Levinger's factor and compare to both the available experimental data and the predictions 
of the local density approximation, based on nuclear matter estimates. The agreement with 
the experiments is excellent, whereas the local density approximation is shown to sizably 
' overestimate the Levinger's factor in the region of the medium nuclei. 
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Within the Levinger's quasideuteron (QD) model [p the nuclear photoabsorption cross section 
■^j- ' aA(E-y), above the giant dipole resonance and below pion threshold, is assumed to be proportional 

Q\ . to the break-up cross section of a deuteron in hadronic matter, ormfiiL): 

o ■ 

£j : ovt(%) = V D crQoiE^) , (1) 

where E^ is the photon energy and Vd is interpreted as the effective number of the nucleon-nucleon 
(NN) pairs of the QD type (see and references therein). Vd is written in the form 
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where A and Z are the mass and atomic numbers of the nucleus and L is the so called Levinger's 
factor. Vd can be calculated for a given nuclear ground state wave function, thus allowing for a 
microscopic interpretation of the phenomenological Levinger's factor. 

The value of L has been extracted from experiments according to the following two models: 
i) the Levinger's model Q, in which ctqd(Ej) is taken as the deuteron cross section damped by 
an exponential function, taking care of Pauli blocking of the final states available to the nucleon 
ejected from the QD, and ii) the Laget's model [H, which associates aQD(E^) with the transition 
amplitudes of virtual (n + /o)-meson exchanges between the two nucleons of the QD pair. 

Both models provide satisfactory fits of photoreaction data in heavy nuclei, but yield different 
values of the Levinger's factor, LL ev (^4) and LL a get(^4)> LL ag et(^.) being ~ 20% larger than LL ev (^4)- 
The effective number of deuteron-like pairs, as well as of three- and four-body structures, in 
spherical nuclei has been investigated within the shell model approach in Refs. 0, ||. In a recent 
paper pj (referred to as I hereafter) we have analyzed the properties of deuteron-like structures in 
infinite symmetric nuclear matter (NM) , described by a hamiltonian containing the realistic Urbana 
v%4 NN potential and the Urbana TNI many-body potential y|. A correlated wave function having 
spin-isospin dependent, central and tensor correlations has been used within the correlated basis 
functions (CBF) theory to compute the QD distribution function in matter and extract the NM 



'Electronic address: Alexei.Illarionov@jinr.ru 



© submitted to 6 th Workshop on "c-m induced Two-Hadron Emission" ,Pavia, 2003. February 8, 2008 - 2 



Levinger's factor at equilibrium density, Ljvm = 11.63, to be compared to the empirical estimate, 
L expt (^ = oo) = 9.26. 

CBF theory has established itself as one of the most effective tools to realistically study, from a 
mi cros copic viewpoint, properties of infinite matter of nucleons ranging from the equation of state 
P. I10j| to the momentum distribution and the one- and two-body Green's functions 1121 1 . In 
the last decade these studies have been successfully extended to deal with finite nuclei 0, 0, ^1 . 

In the paper (referred to as II hereafter) we extended the CBF many body approach, used 
in I for NM, to evaluate ab initio the momentum distribution, Pd(\s.d), and the total number per 
particle, Vd/A, of QD pairs in the doubly closed shell nuclei 16 and 40 Ca, described in the Is 
coupling scheme. From VnjA we then extracted the corresponding Levinger's factors. 

Following the approach developed in I and II, in a ^4-nucleon system the momentum distribution 
of QD pairs can be written as the overlap between the nuclear and deuteron ground state wave 
functions in terms of the natural orbits (NO), <f>a (k) an d their occupation numbers, n a= ( nhn ), 
which are obtained by expanding the one-body density matrix in multipoles. 
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Y/ m (q) being the spherical harmonics. 

In the independent particle model (IPM), an ^4-body wave function ^a(R) is the Slater de- 
terminant of single particle orbitals (fr a (i), which are eigenfunctions of a suitable single particle 
hamiltonian, and n^ M = 1, <p^° = <p a for occupied states, whereas n^ M = for unoccupied 
states. For nuclear matter, the orbitals <f> a (i) are plane waves corresponding to a noninteracting 

Fermi gas of nucleons with momenta |k| < Zip = (6tt 2 pnm /v) 1 ^ , v = 4 and pnm are the NM spin- 
isospin degeneracy and density, respectively. Deviations from IPM provide a measure of correlation 
effects, as they allow higher NO to become populated with n a ^ 0. 

The U(r) = uu(r) — Au(r) and W(r) = wr){r) — Aw(r) functions account for the medium 
correlations effect on the bare components of the DWF, U£){r) and iuo(r). Their explicit expres- 
sions, in terms of the 6 central and non-central spin-isospin depend ent correlation functions, which 
are variationally fixed by minimizing the ground state energy (see |15IJ and references therein), are 
given in I. 



2. RESULTS 



Last generation NN potenti als a re able to fit deuteron properties and the Nijmegen 93 nucleon- 
nucleon scattering phase-shifts up to the pion-production threshold (~ 4000 data points) with 
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a x 2 ~ 1- The Argonne v±s, belonging to this generation, is given by the sum of 14 isosc ala r 
and 4 isovector terms, including charge-symmetry and charge-invariance breaking components |18j |. 
In this work we have used a simpler NN potential, referred to as Argonne v' 8 , obtained from the 
the full Argonne v\s retaining only the first eight operatorial terms, corresponding the central and 
non-central spin-isospin dependent components plus spin-orbit and spin-orbit /isospin components. 
The Argonne v' 8 is constructed in such a way to reproduce the isoscalar part of the full v\% in the 
S, P and 3 .Di waves and the s Di- s S\ coupling. The v' 8 parameterization, while allowing for a 
fully realistic NN interaction, makes the use of modern many-body methods much more practical. 
It has been found that the differences between Argonne v' 8 and the full v\s contribute very little 
to the binding energy of light nuclei and nuclear matter, and can be safely estimated either by 
perturbation theory or from FHNC/SOC calculations. 

It is well known that, to quantitatively describe the properties of nuclei with A > 2, modern 
NN interactions need to be supplemented with three-body forces. The Urbana IX (UIX) model 
provides a very good description of the energies of both the ground and the low-lying excited states 
of light nuclei (A < 8). In the present calculations we use Argonne v' s + UIX interaction, which 
will be referred to as as the AU8' model. This interaction has already been used in the variati ona l 
FHNC/SOC calculations of Ref. as well as i n the quantum Monte Carlo simulations of Ref. |19IJ . 

For the single particle wave functions, (f) a (i), entering the shell model wave function we have 
used the solution of the single particle Schrodinger equation with a Woods-Saxon potential, 

Vws ^ = TT /° bw 1 • ( 8 ) 

1 + exp [(r - Ro)/a \ 

In principle, the parameters of the correlation functions, f p (r), and of the Woods-Saxon po- 
tential may be both fixed by minimizing th e gr ound state energy. This complete minimization 
was performed for the AU8' model in Ref. j!5lj . and provided a binding energy per nucleoli of 

B/A = 5.48 MeV in 16 and B/A = 6.97 MeV in 40 Ca (the experimental values are 7.97 MeV 
in 16 and 8.55 MeV in 40 Ca). These differences are compatible with the results of nuclear mat- 
ter calculations at saturation density, Pnm = 0.16 fm~ 3 , carried out with the same hamiltonian. 
In fact, the FHNC/SOC nuclear matter energy per nucleon is E NM /A = -10.9 MeV [H, to be 
compared to the empirical value of —16 MeV. 

However, the calculated root mean square radii of the two nuclei turned out to be R = 2.83 fm 
in 16 and R = 3.66 fm in 40 Ca, showing a difference of ~ 5% with the experimental values, 
-Rcxpt = 2.73 fm and -R exp t = 3.48 fm, respectively. Moreover, the one-body densities were not in 
close agreement with the experimental ones. In order to take care of this feature of the variational 
approach, a set of single particle wave functions providing an accurate description of the empirical 
densities was chosen, and the energy was then minimized with respect to the correlation functions 
only. The resulting radii were R = 2.67 fm ( 16 0) and R = 3.39 fm ( 40 Ca), with a density description 
very much improved. The energies obtained by this partial minimization procedure were B/A = 
5.41 MeV in 16 and B/A = 6.64 MeV in 40 Ca, largely within the accuracy of the FHNC/SOC 
scheme. Here, we adopt this same wave function, whose parameters are given in Table V of Ref. |15j |. 

The structure of the NO in 16 O and 40 Ca is discussed at length Ref. |14j |. Here we limit 
ourselves to recall some of their main characteristics. The effect of correlations is mostly visible 
in the Is orbital, where the NO are larger than the shell model ones at short distances, resulting 
in stronger localization. The influence on the shape of the other occupied shell model orbitals is 
negligible. The occupation of the NO corresponding to the fully occupied shell model states is 
depleted by 9. 6% in 16 and by 10.5% in 40 Ca, with a maximum depletion of ~ 22% for the 2s 
state in 40 Ca. As a consequence, the lowest mean field unoccupied states become sizably populated 
(ri2s ~ n2p ~ n\d ~ 0.02 in 16 and n% s m 0.05, nip ~ 0.02, n2d ~ 0.03 in Ca). These two effects 
are largely due to the presence of the tensor correlation. 

Fig. n snows the behavior of U(r) and W(r) in 16 0, 40 Ca and nuclear matter, evaluated using 
the hamiltonian AU8'. For comparison, we also show the bare components of the Argonne v' s DWF. 
It appears that the main differences between deuteron and QD occur at r < 2 fm. At small relative 
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distances both U{r) (r < 1 fm) and W(r) (r < 0.5 fm) are slightly suppressed with respect to ud(t) 
and wd{t)- On the contrary, they are appreciably enhanced at larger distances. These effects are 
more visible for the lightest nucleus. 

The differences between nuclear matter and nuclei mostly disappear in the Fourier transforms, 
\U(k)\, \W(k)\, \uo{k)\ and \wD(k)\, whose behavior is displayed in Fig. The nuclear medium 
shifts the second minimum of |u_d(A;)| towards lower values of k, as obtained in I for nuclear matter 
with the Urbana v\4 potential. The Argonne v' 8 \wo(k)\ does not exihibit any diffraction minimum, 
which, however, appears in |W(fc)|. 

The distribution of deuteron pairs with total momentum k^, VdO^-d), resulting from our ap- 
proach is displayed by the solid line in Fig.|3]for 16 and in Fig.^for 40 Ca. The following comments 
are in order: 

(i) NN correlations introduce high momentum components in the distribution. The full VdO^d) 
is strongly enhanced with respect to Vjf M (\tD) at large |koj, and it is correspondingly 
depleted at small |k£>|. The depletion is mostly due to the non-central tensor correlations. 

(ii) The effect of state-dependent correlations is large, as one can see by comparing the full 
VdO<-d) with the Jastrow model "P^ko) (obtained by retaining only the scalar component 
in the two-body correlation operator). 

(iii) The tail of VdO^d) is appreciably different from that of nuclear matter. At |ko| = Akp the 
difference is still a factor ~ 10 for both 16 and 40 Ca. 

Fig. displays the convergence of VdQc-d) in the number of natural orbits included in the sum 
of Eq. ©. Full convergence is reached with the inclusion of orbits up to 5/ for both 16 and 40 Ca. 
The figure shows that, in the case of 40 Ca, the tail of Vd^d) is still ~ 10 times too small if only 
orbitals up to 3d are included. 

The total number of pairs of the QD type in both finite nuclei and nuclear matter, Vd is obtained 
by integration of VdO^-d) over kp: 



where the factor 3 in the r.h.s corresponds to the spin multiplicity of the deuteron, 2Jjj + 1. 

We have repeated the calculations for nuclear matter by using the AU8' interaction. The result 
Vd{NM)/A = 2.707 (the corresponding Fermi gas model result is 3.382) should be compared to 
the value 2.895, obtained in I with the Urbana v\4 two-nucleon plus the Urbana TNI many-body 
forces 8] (which will referred to as the UU14 model). The corresponding numbers for 16 and 40 Ca 
turn out to be much smaller: Vd( 1& 0)/A = 1.090 and Vd( 40 C&)/A = 1.370, respectively. 

The Levinger factor is easily obtained from Vd/A by means of Eq. @. As we are dealing 
with symmetric matter (N = Z = A/2), L(A) = AVd/A. Our estimates L f(3 ( 16 0), L f6 ( 40 Ca) and 
Lf 6 (NM) for 16 0, 40 Ca and nuclear matter, corresponding to the fe correlation model, are reported 
in Figs.|3]and01 These results are not too different from the values obtained within the independent 
particle and Jastrow models. This fact actually implies that the high momentum tail of VdO^-d) 
is not relevant for the calculation of the Levinger factor L. It has to be stressed that the Jastrow 
model turns out to consistently underestimate the Levinger factor. Th e sp atial structure of np pairs 
having the deuteron quantum numbers has been investigated in Ref. 20] in light (A=3,4,6 and 7) 
nuclei and 16 O using a variational Monte Carlo approach and the Argonne v\$ two-nucleon and 
Urbana IX three-nucleon potentials. The estimated Levinger factor for 16 is Lvmc( 16 0) = 4.70, 
comfortably close to our value, L fg ( 16 0) = 4.36. 

Our results for the Levinger factors are summarized in Fig.El where they are also compared with 
the available experimental estimates. The agreement with the photoreaction data of Ahrens et al. 
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22] for the case of 16 and 40 Ca is rather impressive. The "experimental" value, L exp t(oo) = 9.26, 
deduced from the phenomenological formula 

L expt (A) = 13.82-^ , (10) 

reported in Ref. 23], is ~ 15% smaller than our theoretical value. In I, the surface contribution to 
L(j4) has been estimated exploiting the calculated enhancement factor of the electric dipole sum 
rule for finite nuclei, K,, |25j |. obtained using CBF theory and LDA. The enhancement factor is 
related to experimental data on photoreactions through the equation: 

1 + /C cxpt = — / a A {E 1 )dE 1 , (11) 
°o Jo 

where do = 60 [Z(A — Z)/A] MeV mb and is the 7r~meson production threshold. By using the 
same parameterization as in I for the surface term, we get: 

Llda04) = 10.83 - 9.76 A~ 1/3 , (12) 

for the AU8' interaction. Llda(^4) is displayed on Fig. LDA turns out to be not satisfactory 
for medium nuclei, such as lb O and 10 Ca. Fig. GO also report Li, ey (A) and LLagct(-A). as extracted 
[23l . |24] from the available experimental data on photoreactions. The computed Levinger's factors 
are almost ^4-independent for heavy nuclei (A > 100), and result to be ~ 15% larger than LL ev (.A) 
and ~ 25% smaller than L-^ a „ et (A). Such disagreement between theory and experiment is likely to 
be ascr ibed to the sizable tail contributions to the electric dipole sum rule, absent in the definition 
of Eq. (HU). 



3. CONCLUSION 



The Correlated Basis Function theory of the two-body density matrix has been applied to 
microscopically compute the distribution of QD pairs carrying total momentum ko, VdQc-d)-, m - 
doubly closed shell nuclei 16 and 40 Ca and nuclear matter, starting from the realistic Argonne v'g 
plus Urbana IX potential. 

It has been found that NN correlations produce a high momentum tail in VdO^-d) and, corre- 
spondingly a depletion at small k/j for both nuclei and nuclear matter. These effects are mainly 
due to the presence of the state-dependent correlations associated with the tensor component of 
the one pion exchange interaction. Contrary to what happens for the one-nucleon momentum 
distibution, the tail of Vd{^d) sizably differs from that of nuclear matter. 

Summation of VdQc-d) over ko provides the total number Vd of QD pairs, and, consequently, 
allows for an ab initio calculation of the Levinger's factor L(A). The CBF result for nuclear 
matter is significantly reduced with respect to the value obtained in I with the Urbana V14 plus the 
Urbana TNI many-body forces. The corresponding Levinger factors for 16 and Ca are much 
smaller than the nuclear matter value and in very good agreement with the available photoreaction 
data analyzed within the quasideuteron phenomenology. In addition, our results show that LDA 
overestimates L(A) in the region of the light-medium nuclei. 

The Ij(A) resulting from the full calculation are relatively close to the corresponding values 
obtained within the IPM and Jastrow models. Actually, the high momentum tail of VdQc-d) gives 
a small contributions to the Levinger factor. This feature indicates that the approximation used 
in our calculation (which amounts to including only diagrams at the dressed lowest order of the 
FHNC cluster expansion) is fully adequate. However, it should be noticed that the Jastrow model 
underestimates the Levinger factor. 

In addition, the analysis described in this paper shows that when a deuteron is embedded in a 
nucleus, or in nuclear matter at equilibrium density, its wave function gets appreciably modified 
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by the surrounding medium. While in the case of the S'-wave component the difference is mostly 
visible at small relative distance (r < 1 fm), the D-wave component of the QD appears to be 
significantly enhanced, with respect to the deuteron wd{t), over the range r < 2 fm. This effect is 
particularly evident in the lightest nucleus. 
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FIG. 1: Radial components U(r) and W(r) of the AU8' QD wave functions in 16 0, 40 Ca and nuclear matter. 
Upper panel: the solid and dashed lines show the radial dependence of UA(r) for 16 and 40 Ca, respectively. 
The dot-dashed and dotted lines correspond to the nuclear matter Unm{t) and the bare ud{t). Lower panel: 
as in the upper panel for the d-wave components of the QD and deuteron wave functions. 
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L f6 ( 16 0) = 4.3601 ; 
Lj( 16 0) =4.1889 

O) = 4.4585 ! 
L fe (NM) = 10.8272 I 




\k D \/(2k F ) 

FIG. 3: Momentum distribution of QD pairs in 16 as a function of the total momentum |ko| (see Eq. J3J)). 
The solid, dashed and dash-dotted lines are the results obtained within the /g and Jastrow correlation models 
and IPM, respectively. The short-dashed line displays the f@ momentum distribution of the QD in nuclear 

matter at equilibrium density, pnm = 0.16 fm~ 3 . The insert shows a blow up of the region |1cd|/(2A;.f) < 1, 
plotted in linear scale. The Levinger factors, L(-A), for the various calculations are also reported. 

^i 

= 5.4808 
L.j( 40 Caj = 5.4305 - 
Li PM ( 40 Ca) = 5.9947 
L ffi (NM) = 10.8272 I 




|k£,|/(2M 



FIG. 4: As in Fig. for 40 Ca. 
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FIG. 5: Convergence of Vd(^d)/A for 16 and 40 Ca in the number of natural orbits included in the 
summation of Eq. The results have been obtained within the fe correlation model. 
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FIG. 6: Levinger's factor L(A) for 16 0, 40 Ca and nuclear matter (shown by the arrows for the UU14 and 
AU8' forces). The filled circles, the empty circles and the triangles show the Levinger's factors obtained 
within the fe and Jastrow correlation models and the IPM, respectively. The LDA, as discussed in the text, 
is also reported (solid line). The phenomcnological values of LL ev (^4) corresponding to the photoreaction 
data of Lepretre et al. [23 (squares) and Ahrens et al. [23 (crosses and diamonds) are taken from Ref. (23] ■ 
The empirical values of L,Laget (A) , represented by circles in the heavy nuclei region, are from Ref. j23 • 



